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Abstract 



o 
o 

. An instability in the presence of matter in theories of gravity which include a 1 /R 

pLn I correction in the gravitational action has been found by Dolgov and Kawasaki. In the 

l/^ ' present paper this instability is discussed for f{R) gravity in general. We focus on 
the Palatini formalism of the theory and it is shown that no such instability occurs 

! in this version of f{R) gravity. The reasons for the appearance of the instability in 

' the metric but not in the Palatini formalism are fully investigated. 

o ■ 
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O 

"o 

jli ! Modified gravity has been attracting increased interest recently, since it con- 

^1 stitutes an interesting option for explaining the nature of dark energy, which 

according to observation should account for about 70% of the energy density 
^ ■ of the universe and be responsible for its current accelerated expansion [1]. 

What is more, our current candidates for a more fundamental theory which 
will be able to describe gravity in a high energy regime do predict corrections 
in the low-energy effective action for gravity, with respect to the well-known 
Einstein-Hilbert action [2]. 

The simplest form of such a generalized action seems to be 

S = ^l d'x^f{R) + SM{g,u. (1) 



where = SnG, Sm is the matter action, ip collectively denotes the matter 
fields and / is a general function of the Ricci scalar R. The field equations 
of General Relativity (GR) can be derived from this action by varying with 
respect to the metric and applying the stationary action principle, once we 
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assume that f{R) — R. In the most general case one derives the following field 
equation [3]: 

f\R)R^,. - \f{R)9^, - V^VJ'iR) + g^^af = ^T^,, (2) 



where V is the covariant derivative defined with the Levi-Civita connection of 
the metric, T^^ = ^1^' ° = ^^^i^ and A'{x) = dA{x)/dx. 

Starting from the the same action, one can think of treating the metric and the 
connection as completely independent entities. In this case R will be replaced 
everywhere with R = g^^ Rn^jiV), where i?^j,(r) is the Ricci tensor constructed 
in the usual way with the independent connection T^^^^- Assuming that this 
connection is symmetric (torsionless) , the independent variation leads to the 
field equations [4,5] 



nR)Ri,u) - \f{R)9^u = (3) 



'-9 



= 0, (4) 



where V is the covariant derivative defined with the independent connection. 
This constitutes what is often called in the literature, the Palatini formalism. 
Notice that this approach assumes that the matter action does not depend 
on the connection, and is therefore not suitable for all types of matter, such 
as fermions for example, unless one imposes the extra condition that matter 
only couples to the Levi-Civita connection of the metric [5,6]. What is more, it 
would be natural to allow the connection to be non-symmetric if we wanted to 
consider fermions. There exists a natural generalization of the Palatini formal- 
ism, the metric-affine formalism, which does indeed have these characteristics 
and is therefore much more appropriate in the present context [5,7]. In the 
metric-affine formalism there is a matter term on the left hand side of eq. (4) 
coming from the variation of the matter action with respect to the connection. 
However, in the present paper we will restrict ourselves to the Palatini formal- 
ism since this study will focus on instabilities that appear in the weak field 
and therefore we will not need to consider fermions. In this case the Palatini 
and the metric-affine formalisms actually coincide, so our results will apply 
for both cases. 

f{R) gravity, both in the metric and the Palatini formalism, has been widely 
used to address the dark energy problem and it has been shown that models in 
which f{R) includes terms inversely proportional to R (or R) can account for 
the late time accelerated expansion of the universe [8,9]. Note that positive 
powers of the curvature can provide a mechanism for early time infiation 
as well [10]. However such models are not problem free. For example, it is 
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still a matter of investigation whether there can be a smooth passage from 
one cosmological era to the other and whether this can give a cosmological 
evolution including all eras [11]. There is also still debate on which models have 
the correct weak field limit or even which is the correct way to obtain that 
[12]. However, our concern in this paper is to study an instability discovered 
by Dolgov and Kawasaki [13] in the presence of matter for the specific model 
where f{R) = R — fJ'^/R, with /i being a constant [9]. This instability is not 
just a special characteristic of this model but occurs in a more general class 
[14]. 

Let us briefly review the results of [13,14]. By contracting eq. (2) one gets 
3n/'(it:) + f{R)R - 2f{R) = k'T, (5) 



where T — g^'^Tn^,. Following [14] we can write f{R) — R + e(p{R), where 
e is a constant. If we consider a small region in a weak field regime within 
matter we can assume that Qah = Vab + hab and R = —k"^ T + Ri, where riab is 
the Minkowski metric and |i?i/«;^T| <C 1. In this approximation, and to first 
order in Ri, eq. (5) gives 



2^2 /// 

Ri - V^Ri ^ {fiti - VT • VRi) 

where an over-dot denotes differentiation with respect to time, whereas V 
and denote the gradient and Laplacian operators respectively in Euclidean 
three-dimensional space. 

The instability occurs if = f"{R) < and if e is very small, since the 
coefficient of the last term on the left hand side of eq. (6) is the square of an 
effective mass. As already mentioned in [14] it can be considered as an instabil- 
ity in the gravity sector. Because of this, and since it appears in the equations 
governing the dynamics of the curvature scalar, we refer to it as "curvature 
scalar instability". Theories with f"{R) > will be stable irrespective of the 
value of e. However, for several models that lead to the desired cosmological 
dynamics at late times, e is indeed very small and f"{R) is indeed negative. A 
typical example is the model of [9], where (f{R) = —fi^/R, with /x ~ 10~^^eV, 
and the time scale for the instabihty to occur is of the order of 10~^^ s [13]. 

All of the above is with reference to the metric formalism. Let us now consider 
the Palatini formalism. We will argue that such an instability cannot occur in 
this case irrespective of the form of the Lagrangian. Contracting eq. (3) gives 

f'{R)R-2f{R) = ^T. (7) 
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Note that R is not the Ricci scalar of the metric. Eq. (4) can actually take 
the form Vxh^^,^ = if one defines h^^, as = f'{R)g^y, which implies that 
the connections are the Christoffel symbols of this metric [4,5]. We can then 
express R in terms R: 

R^R ^-^^ f(R)^i^f(R) + ^^nf(R). (8) 



Now notice that eq. (7) is an algebraic equation in R for a given f{R) which 
will have solutions of the form R = 0{T), where 6 is some function. We are 
not going to examine cases in which eq. (7) has no solutions or is identically 
satisfied {f{R) oc i?^), since these do not constitute viable choices for a low- 
energy gravitational theory [4,5]. 

We can now write eqs. (8) as 



(9) 

or alternatively R = F{T), where F{T) is a function of T. This clearly demon- 
strates that the Ricci scalar of the metric can be expressed directly as a func- 
tion of the trace of the stress-energy tensor. In fact, eq. (9) is a straightforward 
generalization of the contracted Einstein equation, R = —h?T. From the form 
of eq. (9), it is clear that no instability can occur in this case, since R carries 
no dynamics in eq. (9), unlike eq. (5). 

Let us now analyze where this difference between the two formalisms stems 
from. By generalizing the Lagrangian from R or R one adds inevitably a 
scalar degree of freedom [6]. However this degree of freedom seems to be of a 
different nature in the two versions of the theory. In the metric version it is 
dynamical, and therefore care should be taken to ensure stability, whereas in 
the Palatini version, it is non dynamical. This is related to the fact that the 
Palatini formalism leads to second order equations as opposed to the metric 
formalism which leads to fourth order field equations in the metric (see [5] for 
a discussion), but it can also be seen easily by using the equivalence of f{R) 
gravity and scalar-tensor theory [15,6]. 

Introducing an auxiliary field x, we can write the action for the metric for- 
malism as [15,6] 

s^^Jd'xV^(f{x) + nx)(R-x)) + SM(9,.,^). (10) 



Variation with respect to x leads to the equation x = -R if f"{x) 0) which 
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reproduces action (1). Redefining the field x by $ = f'ix) setting V{^) — 
x($)$ — f{x{^)) the action takes the form 

S = ^j d'x^g{^R-Vm + SM{g,u.^). (11) 

which is the action of a scalar-tensor theory with ujq = 0. In the Palatini 
formalism the action will be formally the same apart from the fact the R will 
become R and $ will become $, but in this case this will not be a scalar-tensor 
theory with a;o = since R is not the Ricci scalar of the metric [6] . However 
if we use eq. (8), and with $ = f'{R): we get 

Spai = 2^ / d^'x,/^ (l>i? + ^d,^d>^^ - l^(l>)) + Sm{9^.u, (12) 

which is indeed a scalar tensor theory, but with cuq — —3/2. 
The field equation of the scalar field in scalar-tensor theory is 

(2u;o + 3) □(/) = K^T + (f)V' - 2V. (13) 

Note that $ and l» are the extra degrees of freedom of f{R) gravity, with 
respect to General Relativity, in the metric and in the Palatini formalism 
respectively. Using eq. (13) it is obvious that they satisfy the field equations 



3n$ + 2V{^) - $y'($) = K^T, (14) 
2V(^) -^V(^)^K^T. (15) 

This demonstrates that $ is indeed dynamical, whereas $ is not, as mentioned 
above. At this point, it is worth mentioning that one should not be misled 
into judging the dynamics of a non-minimally coupled field by the presence 
or absence of a kinetic term in the action. There are no kinetic terms for 
$ in action (11) but it is still dynamical. Exactly, the opposite holds for $. 
The reason for this is that both fields are coupled not only to the metric, but 
also to its derivatives. Therefore, when varying the action with respect to the 
metric and then integrating by parts in order to "free" Sg'^'^, terms including 
derivatives of the scalar field are bound to appear. Therefore, even though 
there are no apparent kinetic terms for $ in the action, there will be kinetic 
terms in the field equations. For $, ujq = —3/2 and this is the remarkable case 
where these kinetic terms exactly cancel out the ones coming from the kinetic 
part of the action. 

It has been shown that the curvature scalar instability discovered by Dolgov 
and Kawasaki for metric f{R) gravity is not present in the Palatini formalism. 
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irrespective of the functional form of /. This is due to the fact that, even 
though the generahzation of the Lagrangian from i? or i? to f{R) or f{R) 
always corresponds to intoducing an additional scalar degree of freedom in the 
theory, this degree of freedom is dynamical in the metric formalism , whereas 
in the Palatini case it is not. This instability does not occur in Palatini f{R) 
gravity but other types of instabilities might well be present. For example, 
judging from the form of eq. (9) it is not difficult to imagine that specific 
forms of / could lead to a blow-up of the scalar-curvature for small density 
perturbations around a stable matter configuration, which are, of course, of a 
different nature. This issue also seems to be directly related to the weak field 
limit of the theory. We hope to address it in future work. 
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